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Abstract. We give a new proof that Rudnick and Sarnak's calculation of the n-correlation 
of zeros of L-functions matches random matrix theory. We do this using the derivation of 
n-correlation for eigenvalues of random unitary matrices from our previous paper. There 
we detailed a method using ratios of characteristic polynomials that is far less elegant than 
standard random matrix techniques, but which has the advantage of translating easily into 
the number theory case by applying the Ratios Conjectures for ratios of L-functions. The 
formulae in our previous paper are unwieldy in comparison with the standard determinental 
random matrix version, but we show here that their form allows immediate identification of 
which terms remain when restrictions are placed on the support of the test function. After 
this it is straightforward to show that the surviving terms match Rudnick and Sarnak's 
result. 



1. Introduction 

Let U (N) denote the group oi N xN unitary matrices and let dU signify the Haar measure 
for this group. For a unitary matrix ?7 we let 6*1, . . . , denote the N eigenangles, i.e e*^^ 
are the eigenvalues. A beautiful theorem (see [Meh]) in random matrix theory gives the 
n-correlation of these eigenvalues: 

Theorem 1. Let F : [0,27r]" C be an integrable function of n-variables. Then 

I ^(%, • • • , OjJdU = -J— I F{ai, . . . , a„) det S^iak - aj) da 

Ju{N),. {^V J[0M" """" 



1 • • • WUji, 



dOr 



indicates that the sum is over n-tuples of distinct indices (ji, . . . , j„) with 1 < 
ji < n, i.e. ji 7^ jk if i ^ k, and where 

sin — 

(1) Sr^ia) - ' 



sin I 



In our previous article [CoSn] we give a different derivation for the n-correlation for the 
eigenvalues of matrices in U (N) that is far less elegant, but which can be copied step by step 
to obtain the n-correlation of zeros of L-functions subject to the Ratios Conjecture. 
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It is of interest to ask about 




and the scaled limit of this quantity as — )■ oo in the situation where / is a translation 
invariant function that has a Fourier transform with limited support. Such a quantity arises 
in number theory, for example in the calculation by Rudnick and Sarnak [RuSa] of the n- 
correlation of the zeros of the Riemann zeta-function. The inherent difficulty of dealing with 
off-diagonal terms in the mean-value theorems which arise necessitate the limited support 
of /. These off-diagonal terms lead to the consideration of the distribution of prime pairs in 
short intervals. See the work of Bogomolny and Keating [BoKe] for a beautiful demonstration 
that the combinatorics of prime number sums arising from the Hardy-Littlewood prime pair 
conjectures lead to the determinantal formula of random matrix theory. 

In [RuSa], after deriving an expression for the n-corrclation in the situation that their 
test function has a Fourier transform with support limited to a proper subset of [—1,1], 
the authors are faced with the non-trivial task of verifying that their answer agrees with 
theorems from random matrix theory. A similar problem occurs in two works on the n-level 
density of zeros of quadratic L-functions, first considered by Mike Rubinstein [Rub], and 
subsequently by Peng Gao in his thesis. In the former case Rubinstein finds, as do Rudnick 
and Sarnak, an adhoc method to verify the consistency of the number theory and the random 
matrix theory calculations. Peng Gao is unable to make this verification in his situation. A 
subsequent paper [ERR] rectifies this situation and completes the verification in Gao's case 
by a very clever appeal to zeta-functions over function fields. 

In this paper, we show that the formula for n-correlation derived in [CoSn] allows for 
immediate identification of the terms which survive a restriction on the support of the test 
function, and so makes the verification of Rudnick and Sarnak's formula in the unitary case 
straightforward, eg. for the n-correlation of the zeros of the Riemann zeta-function. (In 
[RuSa] more general L-functions are considered; here we focus on the case of the Riemann 
zeta-function, but higher degree L-functions don't present any extra difficulty.) Our proof is 
natural in that it explains the situation for a test function whose Fourier transform has any 
range of support. The extension of our method to orthogonal and symplectic ensembles will 
be part of Amy Mason's thesis. 



2. Statement of results 

For simplicity we state the results of Rudnick and Sarnak for the n-correlation of the 
zeros 1/2-1- i'fj of the Riemann zeta-function, though there is no difficulty working in their 
full generality. Before doing so, we introduce their vector notation as a convenient way to 
express the combinatorial sum that arises in their work. Let 

I — ei — ej 

^ ^ \ ei = (0, . . . , 1, . . . , 0) the ith standard basis vector 
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Theorem 2. [RuSa] Theorem 3.1. Let hj, for I < j < n, be rapidly decaying functions with 

(4) hjix) = [ ^,(i)e-* dt 
where Qj is smooth and compactly supported. Suppose that 

(5) f{xi,...,Xn) ^ $(^1, . . H hCn)e(-Xi^i XnCn)dCl ■ ■ ■ d^n 

where $ is smooth, even, and compactly supported in such a way that 

(6) $(ei,...,U = 

whenever 

(7) |ei| + --- + |en|>2-e 
for some fixed e > 0. Let C — logT. Then 



7l,-,7n 

[n/2] 

= K 



27r 



TC ( f \ 

(h) 2ir( *(0) + J2 J2 l^il ■ ■ ■ + • • • + ^r-ei(.)j(.)) dv\ + 0{T) 

V r=l i(t)<j(t) ^ 

t<r 

where the sum is over all disjoint pairs of indices i{t) < j{t) in {1,2, ... ,n} when t <r and 
where 

(8) «;(h) = / hi{u) . . .hn{u) du. 

The main resuh of this paper is the random matrix analogue of this theorem. For the 
purposes of making a close connection with number theory it is convenient to let the eige- 
nangles "wrap around." Thus, for a,n N x N unitary matrix X we will let its eigenangles 
be 

(9) ■■■<e-R<0 -R+i ^ • ■ ■ ^ Oq < 9i < ■ ■ ■ < 9r < . . . 
where 

(10) Or+kN ^Or + 27rk. 

Now we have infinitely many eigenangles just as for the Riemann zeta- function we have 
infinitely many zeros. We are interested in the sum 



Jl,-,Jn 

where each jk now runs over all integers. This sum exactly parallels the above sum over 
zeta-zeros. 
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Theorem 3. With the same conditions on h and f, we have 



(12) E E / e.....e..$(Eefc,e,^,,^^Jde + 0(7V)- 

K+L+M = {l.. .,u} aeSx ^''j-^^ \j=l / 

\K\ = \L\ 

This theorem is shghtly more general than the theorem above of Rudnick and Sarnak and 
has a smaller error term; but the main terms exactly match when $ is assumed to be even 
as in [RuSa] and N takes the place of logT". As a consequence, we have proven that the 
result of Rudnick and Sarnak agrees with random matrix theory without going through the 
complex combinatorial considerations that they undergo in the second half of their paper. 

In general, one could ask to compare number theory and random matrix theory to see 
what happens if one works with a set of test functions Tg which are defined as above except 
that the support of $ is restricted so that 

(13) $(ei,...,en)=0 

whenever 

(14) \^i\ + --- + \U\>'2q-e. 

We have developed an approach in [CoSn] that treats number theory and random matrix 
theory in parallel. Prom this work it is a simple matter to read off the results when dealing 
with test functions from Tq. This result is described below. 

3. Eigenvalue correlations 

Before stating the theorem of [CoSn] we describe some notation. We let 

1 



(15) z(x)^ _^ 

1 — e ^ 

In our formulas for averages of characteristic polynomials the function z{x) plays the role 
for random matrix theory that C(l + 2;) plays in the theory of moments of the Riemann 
zeta-function. 

Given finite sets A and B we will have a sum over subsets S G A and T G B with 
IS"! = |T|. We let S = A — S and T = B — T. Wc will let a denote a generic member of S 
and /3 denote a generic member of T; we will use a and /3 for generic members of A and B 
or of S and T, according to the context. Also S~ — {—a : a E S}, and similarly for T~. We 
let 

(16) Z{A,B):^llz{a + (3). 
A simple modification of Theorem 4 of [CoSn] yields 
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Theorem 4. Let S > and let J^^^ denote an integration along the vertical path from c — ioo 
to c + ioo. Suppose that F{xi, . . . is a holomorphic function which decays rapidly in 
each variable in vertical strips. Then, for any 5 > 0, 



I 



J] F{e,„...,e,jdx 

(17) =— V F (_iyL| + |M|^|M| 

^ ' K+L+A^={l,...,n} 

X / / / J*{zk\ -ZL)F{izi, izn) dzi... dz^ 

J{-8)\L\ J(o)\M\ 

where zk = {zk '■ k G K}, —zl = {—zg : i E L} and J^^^yKi /(o)i*f| t^^ans that we are 

integrating all the variables in zk along the (5) path, all of the variables in zl along the (—5) 
path and all of the variables in zm along the (0) path; and 

J*{A,B) : = 

(18) V e-^(^<^^5"+^/^e^^) Z{S,T)Z{S~,T-) f]H.T(U) 

SCA,TCB ^ ' / V ' / (A-S) + (B-T) o=l 

\S\ = \T\ =Ui+- + Uq 

\Ugl<2 

where 

f E..5T(«-«)-E^eT7(« + /3) ^fW = {a}cA-S 



(19) Hs,t(W) - <^ 



{a + /3) tfW = {a, with ^^^ij 

otherwise. 



Also, Z{A,B) = Yl^eA z{a + (3), with the dagger imposing the additional restriction that a 
factor z{x) is omitted if its argument is zero. 

We remark that, assuming the ratios conjecture, see [CoSn], there is a structurally identical 
formula giving the n-correlation of the zeros of the Riemann zeta-function. 
We want to apply the above but with 

(20) .(..,...,.„).,(^..,„i^)m)....„(| 

with f GTq. 

In the right-hand side of the formula we replace / by its Fourier transform so that we can 
better see what the implications of limited support are. Thus, we write 

/H^^'---'^^^ = / *(6,---,en)(^(6 + --- + Ue ^- 



27r 27r / ./r™ V 27r 27r 



6 J.B. CONREY AND N.C. SNAITH 

Observe that 

(21) e ^Nzn^n ^ ^ ^Nz,i^+-+Nz„^„_ 

For each i, \ReZi\ < S. Suppose that $(^i, . . . , ^n) = if |^i| + ■ ■ ■ + > 2q — e for some 
e > 0. Then |$| 7^ imphes that 

(22) |gA^-2l5i + ---+A^2n5n I ^^N5{2q-e) 

We compare this exponential with the exponentials which appears in the factor J*{zk', —zl)'- 



(23) E 



\S\ = \T\ 



Notice that the real parts of all of the a E zk and all of the /3 G —zl are equal to 6. If 
15*1 = |T| > q, then 



(24) 



Thus, the product of these two factors is < e"^^"^. We can move the paths of integration 
(—5) and (S) away from the imaginary axis. The integrand tends to zero uniformly on the 
vertical paths as 5 — )■ 00. Thus, we see that these terms are 0. 

We let Jg{A; B) be defined as J* [A] B) but with the subsets S and T in the defining sum 
having size smaller than g, i.e. 

(25) rM-,B)= Yl ■■■■ 

SCA,TCB 
\S\ = \T\<g 

Then, if the total support of $ is limited to any number smaller than 2q, then Theorem H] 
holds with all the J* replaced by J*. 

Theorem 5. Let 6 > and suppose that F{xi, . . . , Xn) satisfies (2U\) . Then, 

(26) =^-V y (_i)|i'l+lM|;vl^| 

^ ' if + L + M = {l,...,n} 

|A'| = |L| 

/ JI{zk] -ZL)F{iZi, izn) dzi... dzn 
(<5)l-»l J{-^5)W J{0)\M\ 
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4. The special case g = 1 



If g = 1, as in the theorem of Rudnick and Sarnak, then the sets 5* and T in the sum 
defining J* {A; B) are both empty. We have 

Q 

(27) Jl{A-B)= \{HUU,). 

A+B= q=l 
Ui + - + Uq ^ 

\Uq\<2 

We have 

(28) H,,iU) = I + if = ^} -ith 



otherwise. 
This gives, when A = {ai, . . . , ax}, B = {f3i, . . . , /3k}, 



(29) JKA; B)=J2Il (7) ("^ + M- 
Thus, we have 

Theorem 6. Let S > and suppose that F{xi, . . . ,x„) satisfies (2U\} with q = 1. Then, 

(30) =— ^ y (_l)\L\+\M\^\M\ 

^ ' X+L+M={l,...,n} 



where 

(31) K = {h,...,kK} and L = {h, . . . , El}- 

5. Comparison with Rudnick-Sarnak 

We now put the formula from the last section into the form of Theorem 3.1 of [RuSa]. 
For ease of subscripting, consider one particular term, denoted (with an abuse of notation 
that hopefully won't cause confusion) K = {1,2, . . . , K}, L = {K + 1, . . . , 2K} and a is the 
identity permutation. Let 



/ / / TT ( ~ ) (^^ ~ ZK+k)F{iZi, izn) dzi... dzn. 
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We write 

x^(^i + ■ ■ ■ + in)e I 1 rf^l . . . rf^„. 

We integrate the Zr variables with 2K < r <2n using equation (jlj) to get 



(0) 



n-2K 



e ^ ^ J /^2x+i J • • • ^« J dz,K^^ ...dz^ 

= n / f f ) ^ f ^) = n (-2-T,,(iVT0)). 

,i;^+ii(o) v^/ V 27r ; ^.^A^^ 

This gives 

1= + + rr {-2mTg^{NTQ) / TT 

-iNzkCk\ , f izk\ f iNzk+K^K+k 
hk [ — \ e 



2tt J \T J \ 271 

X hx+k (^—^-^^ (^-^^ i^k - ZK+k)^ dzi-.. dz2Kd^. 

Now we move the paths of integration of the with 1 < k < K. If .^^ > we move the path 
to the left across the double pole at z^+k] if ^fc < then we move the path far to the right. 
We use the fact that 

fl{Zk)f2{zK+k) ell \f ( \ 

~l ^ = Jl{ZK+k)j2{ZK+k). 

Zk=ZK+k [Zk-ZK+k) 

In this way we get 

/n 
<f(6,..-,U5(ei + --- + en) n i-^rrzTg,{NTQ) 

tj>0,3<K j=2K + l 

K 

X 



jj ( ^^^^ ( + /Hp) A,,,, (!fp) ,,,,,, ) 



I k=l 

We compute, for example 
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Thus, 



/n 
6---e/^^(6,---,en)5(6 + --- + en) n (-27rir^,(A^Te,)) 

ij>0,3<K j=2K+l 

X n f-27riT / 9k{uk)9K+k{-Uk + TN{^k + U+k)) duk] d^{l + 0{^)). 
k=i V / 

We make the changes of variables 

yj = NT^j for 27f + 1 < J < n 

Vk = Uk for 1 < k < K 

VK+k = -Uk + NT{^k + iK+k) iox l<k<K. 
The last substitution implies that 

(32) e,,, = _& + (»£±i±ii5). 

Also, the condition YTj=i ij — imphed by the delta-function translates to Vj ~ 0. We 
have 



{NTy-^-^ 



-K 



10 n 



^k>0,l<k<K 



. . , VK+l + yi ^ , {.yiK + VK) yn-2K Vr. 

X*l?l, ■ ■ ■ , -Ql H 7777^ , ■ ■ ■ , -QK + 



TVT ' ' NT ' NT ' 'NT' 

X (1 + 0(^))(i6 ...d^K dyi--- dyn- 
Employing the Taylor expansion of just as in [RuSa], we have 
I 



N 



2K-n 



{2mn-ir-^ llgj{yj)dyl Ci-.-^ 

E!/,=0 i=l «fe>0,l<fc<K 



NT 

x$(6, ...,^K, -a, • • • , -Ci^, 0, . . . , 0) d6 . . . + ^(l/T)). 



Since 



n—l -, « n— 1 



„ n^j(^.?)^^ ^ / Yl9j{yj)7j- hn{t)exp{-it{J2yj) dt dy 
(34) ^ _L/J.^(,)...^. 



we have 



NT V / V / 27r / 

?fc>0,l<*:<if 



x$(ei, . . . , Ck, -6, ■ ■ ■ , -Cx, 0, . . . , 0) ■ ■ ■ + o(i/r)). 
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More generally, for K = {ki, . . . , k^K\}, L = {ii, . . . , i\K\} and a G S\k\, 
I{K,L,a): = [ [ f f\ ( —) {z,,^ - )F{izi, . . . ,iZn) dz 

i(5)|/v| J^_s)\L\ J(o)|M| yz J 

We insert this into Theorem [6] and have 



1^ y (_l)l^l+|MliVl^^liV2l^l-"(27rO"(-l)"-l^l 

rrrf \ TL ^ ^ 



^ ' K+L+A/={l,...,n} 



\K\ = \L\ 

K(h)NT 

X 



27r 

Since \M\ = n — 2\K\ the right-hand side simplifies to 



NT r \ 

if+L+j\/={i,...,n} o-eSiKi \ j=l / 



\K\ = \L 

This proves Theorem [31 
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